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SUMMARY 


A numerical method is presented for the stress analysis of 
stiffened-shell structures of arbitrary cross section under nonuniform 
temperature distributions. The method is based on a previously pub- 
lished procedure that is extended to include temperature effects and 
multicell construction. The application of the method to practical 
problems is discussed, and an illustrative analysis is presented of a 
two-cell box beam under the combined action of vertical loads and a 
nonuniform temperature distribution. 


INTRODUCTION 


The effects of nonuniform temperature distributions, such as those 
produced by aerodynamic heating, are becoming of greater concern in the 
design of modern high-speed aircraft. The structural effects of tem- 
perature changes and the results of some analyses of a simplified struc 
ture under nonuniform distributions of temperature have been discussed 
in reference 1. The analytical methods considered in reference 1 were 
found, however, to yield inaccurate values for the secondary stresses 
in complicated structures, and in such cases some type of numerical 
approach is desirable. Numerical methods, however, usually require 
extensive and tedious calculations and they should be used only when 
satisfactory results cannot be obtained from a simplified analysis. 

Several numerical methods of stress analysis have been presented 
in the literature, but none contains provisions for temperature changes 
In the present paper, one such method, the numerical procedure of refer 
ence 2, has been extended to include the effects of a nonuniform dis- 
tribution of temperature. In addition, the equations developed permit 
the analysis of a stiffened-shell structure of arbitrary cross section 
with any number of internal cells. The application of the method is 
discussed and illustrated by analysis of a two-cell box beam under the 
combined action of vertical loads and a nonuniform temperature 
distribution. 
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DESCRIPTION OF THE NUMERICAL METHOD 
Basic Theory 


The structure analyzed is an idealized representation of a multi- 
cell stiffened-shell structure (see fig. 1) and has the following 
characteristics: 

(1) The basic unit is a rectangular panel bounded on two parallel 
sides by extensionally flexible stringers and on the other two sides by 
rigid bulkheads. 

(2) The panels consist of sheet material and are assumed to carry 
shear stress only. The shear stress is constant within a given panel. 

(3) The stringers run parallel to the direction of the primary 
stresses and carry axial load only. 

(4) The bulkheads lie perpendicular to the stringers and are rigid 
in their own plane but' offer no resistance to warping out of their plane. 

(5) The structure is loaded only at the bulkheads. 

(6) Material properties, cross-sectional dimensions, and tempera- 
ture distribution do not vary along the length of a given bay. 

With these assumptions about the basic elements of the structure, 
any type of stiffened shell can be analyzed, provided taper is excluded. 
The state of stress in such a structure can then be defined by suitable 
stress-strain relations and two types of displacements: 

(1) Stringer displacements, which are displacements, at the end of 

a bay, of each flexible stringer in a direction parallel to the stringer. 

(2) Bay displacements, which are translations and rotations of the 
plane of each cross section defined by the rigid bulkheads. 

Once the stress-strain relations are established for the components 
of the idealized structure , equations of equilibrium can be used to 
obtain relationships between the displacements. The equilibrium equa- 
tion for the forces on an individual stringer yields an expression for 
the stringer displacement of any panel point in terms of the surrounding 
stringer displacements and the displacements of the two adjacent bays. 
From this general expression, equations equal in number to the unknown 
stringer displacements are obtained. The additional equations required 
for the determination of the bay displacements are obtained from the 
equations of equilibrium of the shear forces on the cross sections. 
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These equations then completely define the displacements of the struc- 
ture. In most cases the number of equations is so large that a direct 
solution would be impractical and it has been found expedient to solve 
them by the recommended iteration procedure described in the next sec- 
tion. The required equations are derived in detail in the appendix. 


Solution of Equations by Iteration 

Matrix iteration often provides the easiest and quickest solution 
to the equations, and the procedure recommended is as follows: 

The equations to be solved can be written in matrix notation as 


For purposes of iteration, these equations are rearranged to give 

{•}-[#}•{•} 

where 

DO - DO - M 



square matrix of coefficients of general equations with diagonal 
terms reduced to unity 



H 

B 


unit matrix 

column matrix of stringer and bay displacements 

column matrix of constant terms in general equation; these terms 
arise from applied load and thermal expansion 


Initial approximate values of stringer and bay displacements -jd 0 f 

are then selected. These values may be determined in any convenient 
manner; however, subsequent operations can be simplified, as explained 
in the appendix, if these values are chosen to correspond to elementary, 
theory. Next, the initial displacement values are substituted into the 


right-hand side of equation (2) to obtain a second approximation 



V 
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to the displacements 



and the differences between the second and initial approximate displace- 
ment values are computed from the equation 




The iteration 
differences. The 


process is then begun by using these displacement 
nth difference is defined as 


|Ad n j = jd n j - jdoj (5) 

and it can be easily verified that the use of these differences leads to 
the following matrix equation: 


|^nj = [ C ]| M n-l| + j^l 


( 6 ) 


The iteration process consists of a series of solutions of equation (6), 
each successive solution yielding a better approximation to the displace- 
ment differences than the previous one. The process is continued until 
successive solutions of equation (6) yield the .same result, that is, 
until 



( 7 ) 


The final displacements are then determined from the final differences 
by using equation ,(j?) and the initial values. 


When equation (6) is being iterated, improved values should be used 
as soon as they are obtained) that is, each individual difference Ad n 


should be substituted into the -} Ad n _q 


r matrix immediately after calcu- 


lation rather than at the end of the cycle. In this manner, each new 
value determined receives the benefit of all previous work and convergence 
is speeded. 


/ 
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The iteration of differences reduces the work required to obtain a 
solution because smaller numbers are involved. However, it is essential 
that no errors be made in the determination of the first differ- 
ences -jAd^ since a single significant error will render the whole 
solution useless. 


Convergence of the Iteration Process 

Li order to obtain more rapid convergence, of the iteration process, 
bay displacements and loads are referred to the principal shear axes of 
each bay. The use of these 'axes greatly simplifies the equations for 
bay displacements by making each bay displacement independent of all 
other bay displacements and thus a function of the stringer displace- 
ments only. In addition, a special correction cycle is periodically 
introduced to bring the stringer forces on each cross section into equi- 
librium with the applied loads. Mathematically, the correction cycle is 
a special cycle that uses a certain combination of the basic equations. 
Its success in the particular case of the numerical method of stress 
analysis is a result of its physical significance, and in that respect 
it is similar to Southwell ’ s "group relaxations" (reference 3). 

The optimum frequency of application of the correction cycle depends 
largely on the characteristics of each individual problem and must be 
determined on a basis of experience with the method. If this frequency 
cannot be determined from previous experience, it can be approximated 
satisfactorily by one that permits the disturbances to spread their 
significant effect over the structure between correction cycles. 

The application of the correction cycle begins at a station where 
the displacements are known and then proceeds outboard. The corrections 
required to bring the first bay into equilibrium are determined and the 
stringer displacements at its outboard end are changed accordingly before 
the corrections required by the second bay are calculated. 


Effect of Introducing Nonuniform 
Temperature Distributions 

The preceding method is applicable to any type of stress problem. 
Nonuniform temperature distributions do not affect the general procedure 
but merely change the details. These effects are of two types: A 

change in the effective structure due to changes in elastic properties 
of the material with temperature and thermal stresses resulting from 
restrained thermal expansion. The changes in elastic properties are 



6 


NACA TN 22kl 


easily handled if the moduli are treated as variables during the deriva- 
tion of the equations. Their effect is analogous to that of variations 
in stringer area and panel thickness. The presence of thermal expan- 
sion requires modification of the stress -strain relationships for the 
stringers but does not affect those for the panels. The equations for 
stringer displacements contain thermal -expansion terms that are analogous 
to the applied-load terms. Bay-displacement equations are unaffected by 
thermal expansion, but thermal -expans ion terms appear in the equations 
used for the correction Cycle. If a difference solution is iterated, 
the elementary solution should include the distributions of thermal 
strain associated with the primary thermal stresses, which may be 
obtained from the equations derived in the appendix. 


DISCUSSION OF THE NUMERICAL METHOD 


The application of a method, such as that just described, always 
poses a number of questions; for example, what are some. of the limita- 
tions of the method, would it be advantageous to use some other method 
of analysis, and how should the structure be idealized. Some of the 
factors requiring consideration, other than those mentioned in the pre- 
vious section, are therefore now discussed. 


Validity of Basic Assumptions 

The assumptions upon which the method is based are commonly 
accepted in the analysis of stiffened shells. Comparison of theoretical 
and experimental results has established the fact that these assumptions 
will yield good results inmost cases. Two important assumptions, that, 
the bulkheads are rigid in their own plane and that the shear stress is 
constant in a given panel, may, however, introduce significant errors 
into the analysis in some cases. These assumptions are therefore 
examined in detail. 

The assumption of rigid bulkheads is satisfactory as long as the 
primary stresses run perpendicular to the bulkheads, but, as demonstrated 
in reference 1, this assumption may not be good when dealing with prob- 
lems involving thermal stress. Large temperature gradients along the 
length of the structure or across the depth of a bulkhead distort the 
real bulkhead and make the assumption of rigidity inapplicable . In many 
cases, however, these effects are small and the assumption yields 
satisfactory results. 

The numerical method could be extended to include the effects of 
bulkhead flexibility. Such an analysis, however, is very cumbersome and 
tedious and if the equations are solved by iteration the process is 
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°f ten very slowly convergent. Therefore, these extensions are not 
considered herein. 

The assumption of constant shear stress in a given panel simplifies 
the development of the equations, and it yields satisfactory results if 
the bulkheads ai*e . reasonably close together. Cases arise, however, in 
which the assumption will lead to unreasonable results because the 
assumed constant shear stress is a poor approximation to a shear stress 
which should be changing rapidly in the spanwise direction. This situa- 
tion is usually accompanied by slow convergence of the iteration process. 
This difficulty, however, can be minimized by reducing the bulkhead 
spacing of the idealized structure since it occurs only when the total 
shear stiffness of the panels joined to a stringer exceeds the exten- 
sions! stiffness of that stringer. 


Idealization of an Actual Structure 

The idealization process described in reference 2 is straightforward. 
However, it provides an opportunity for the stress analyst to exercise 
his engineering judgment and thus simplify the analysis. By restricting 
the analysis to only a part of the structure or by using a comparatively 
simple idealized structure, the time required for the analysis can be 
substantially reduced. Such simplifications, however, can reduce the 
value of the results, and a compromise between speed and exactness is 
required. 

The number and location of the idealized stringers completely 
define the stress-distribution shapes obtainable from the analysis. 

(For example, in an idealized shell of n stringers, there are n pos- 
sible types of independent normal-stress distributions, three of which 
can be determined from elementary theory, the remaining n - 3 being 
statically indeterminate.) Stringer location is thus an important part 
of the idealization process and in conventional problems the locations 
should be selected after consideration of the characteristics of the 
actual structure, the nature of the expected results, and the time avail- 
able for the analysis. When nonuniform temperature distributions are 
involved, the shape of the temperature distribution should also be con- 
sidered because the thermal-stress and temperature distributions will 
have similar shapes and the analysis will yield good results only if the 
idealized structure permits a stress distribution of that shape. 

The bulkhead spacing usually is the same in the idealized and 
actual structures, but the idealized spacing should never be so large 
that trouble is caused by the assumption of constant shear stress. A 
proper bulkhead spacing is one for which the extensional stiffness of 
each stringer element is greater than the sum of the shear stiffnesses 
of the adjacent panels so that no negative terms appear on the right-hand 
side of the equations for the stringer displacements. 
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Calculating Procedures 

All the calculations required by the numerical procedure (deter- 
mining the coefficients of the equations, solving by iteration, and 
computing the stresses) are routine and involve only simple arithmetic. 
The calculations can be easily arranged in tabular form so that the bulk 
of the work can be done by modern automatic computing machinery or by a 
computer who does not need to have a knowledge of the structural 
theories involved. 

In any problem that involves extensive numerical work, errors are 
very apt to occur. One of the advantages of the numerical method 
described herein is that a number of checking procedures can be devised 
to check the various steps in the calculations. No attempt is made to 
describe the many possible checks; a few, however, have been indicated 
in the illustrative example. 

Solution of the equations by simple iteration also possesses 
another advantage with regard to errors. Values obtained from succes- 
sive cycles of iteration show trends that can be observed by an experi- 
enced computer, and errors can be detected by their effects on these 
trends. Errors that do appear during the iteration process eventually 
work themselves out but may adversely affect the rate of convergence. 

APPLICATION OF THE NUMERICAL METHOD 
Description of the Problem 


The application of the numerical method is illustrated by an 
analysis of the idealized two-cell box beam shown in figure 2. The 
cross section is symmetrical about the horizontal center line and the 
beam is untapered; however, the stringer areas and sheet thicknesses 
vary from bay to bay. 

■ The box beam is loaded by four concentrated vertical loads applied 
at the bulkhead stations along the inner web; in addition, it is sub- 
jected to the arbitrarily selected distribution of temperature increase 
shown in figure 3- The temperature is highest at the tip and along the 
front web and decreases in both spanwise and chordwise directions, but 
it is constant across the depth of the beam. The beam is assumed to be 
constructed of 75S-T6 aluminum alloy which has the variation of elastic 
properties with temperature increase shown in figure 4. These data are 
the same as those used in reference 1. 

It is assumed that no thermal stresses were present at 60° F. 

Since the method of analysis involves the assumption that no changes in 
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temperature distribution occur over each element, the temperature used 
in the calculations was the temperature at the center of the element 
concerned. 


Details of the Analysis 

Since the structure and the temperature distribution are symmetri- 
cal about the horizontal center line, the analysis can be restricted to 
one cover. Two solutions are required, however, to determine the total 
stress since the thermal -stress system is symmetrical about the hori- 
zontal center line but the load-stress system is anti symmetrical. In 
this way, the analysis requires the solution of two sets of equations 
(one of 20, the other of 24) which can be solved more easily than the 
set of 44 needed for a single analysis of the complete box. 

The computations required are given in tabular form with most 
tables containing two parts, one related to the load stresses and the 
other related to the thermal stresses. The final solution is obtained 
by the superposition of these two solutions. The rectangular cross 
section and its symmetry permit several simplifications of the general 
equations. In each case the equations used are listed. The notation 
is described in the appendix. Methods used to check the calculations 
are also given in the tables. The checking methods used were determined 
from mathematical relationships existing between the coefficients of 

the equations and from equilibrium of forces. 

/ 

Tables I and II present the physical characteristics and stiffness 
•parameters of the individual stringers and panels. Table III gives the 
location of the principal shear axes of each bay and .the coefficients of 
the bay-displacement equations. The location of the principal inertia 
axes of each bay, the coefficients used in the correction cycle, and 
the initial stringer displacements, are given in table IV. The coeffi- 
cients of the stringer displacement equations are tabulated in table V. 

Table VI contains the matrices used for the iteration. The rows 

and columns have been interchanged in order that the matrix multiplica- 
tions required will consist of the cumulative multiplication of the 
adjacent numbers in two columns. Table VII is a similar arrangement of 
data required for the correction cycle and also lists each correction 
determined. The displacements Obtained from each cycle of iteration 
are given in table VIII and the correction cycles are indicated. 

Table IX contains the calculation of each type of stress and the super- 
position required to obtain the total stresses. 

The numerical calculations in this example were done by a computer 
who had previous experience with the method. The following times were 
required: . 
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Setting up the equations (table I to table VII ) 3 days 

Solving the equations (table VIII) .4 days 

Computing stresses (table IX) 1 day 


In this example, the displacements were computed to six decimal 
places (five or six significant figures) in order that the stress would 
be accurate to 1 psi and thus would provide accurate equilibrium checks. 
Most practical problems will not require such numerical accuracy and a 
smaller number of decimal places should be used -in order to speed the 
solution. It is estimated that the time required to solve this example 
could have been reduced by one-half if the number of dec imal places had 
been reduced from six to four. This reduction would have given stresses 
accurate to 100 psi or about 1 percent of the max imum stress. 


Results of the Calculation 

The results of the calculation are shown graphically in figure 5 
by spanwise and chordwise plots of the stringer stresses in the top and 
bottom covers and a spanwise plot of shear stresses in the webs. The 
spanwise plots have a jagged appearance because stringer areas and sheet 
thicknesses are assumed constant in each bay with an abrupt change at 
the bulkheads. The dashed lines in the plots of stringer stresses are 
the values obtained from an elementary analysis. 


CONCLUSIONS 


A numerical method for the stress analysis of stiffened-shell 
structures under nonuniform temperature distributions has been presented. 
The method is not applicable, to the solution of all structural problems 
involving temperature effects because it requires extensive and tedious 
calculations and because the basic assumptions of bulkheads rigid in 
their own plane and constant shear stress in a given panel occasionally 
lead to unsatisfactory results. It is, however, a powerful tool for the 
solution of many structural problems because: 

(1) It is a means for accurately determining all types of secondary 
stresses in complicated structures that cannot be satisfactorily 
analyzed by simplified methods. 

(2) It is sufficiently flexible to cope with a wide variety of 
structural problems involving nonuniform temperature distributions. 
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(3) It involves only simple arithmetic that can be handled by 
automatic computing machinery. • . 

Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., September 12, 1950 
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APPENDIX 

DERIVATION OF GENERAL EQUATIONS 

The general equations required for the numerical analysis of a 
stiffened shell of arbitrary cross section with any number of internal 
cells and under a nonuniform temperature distribution are developed. 

The basic assumptions and a general description of the method have been 
given previously and are not repeated. 

Symbols 

A cross-sectional area of stringer, square inches 

b width of panel on k grid line, inches 

E modulus of elasticity, psi 

F applied force, pounds 

G modulus of rigidity, psi 

h width of panel on ^ grid line, inches 

I moment of inertia, inches^ 

J shear stiffness parameters 

l,m coordinates of a special set of axes 
L length of bay, inches 

M . ' applied moment, inch -pounds 

P axial load in stringer, positive for tensile load, pounds 

Q area moment, inches^ 

r normal distance to panel on k grid line, positive in positive 

z-direction, inches 

T temperature increment, measured from temperature of zero thermal 

. stress which is 60° F in the example presented, degrees Fahre nhe it 

t panel thickness, inches 
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u, v,w displacements in x-, y-, and z-directions, respectively, inches 

x,y,z rectangular coordinate* axes 

a coefficient of thermal expansion, inches per inch per degree 

Fahrenheit 

P angular rotation used in correction cycle, radians 

y shear strain, radians 

5, A increment 

€ normal strain, inches per inch 

9 angular rotation about x-axis, radians 

X. rotation of special set of axes, degrees or radians 

p normal distance to panel on grid line, positive in positive 

y-direction, inches 

a normal stress, positive for tensile stress, psi 

r shear stress, positive in direction of associated coordinate 

axis when tensile stress on cross section is in positive 
x-direction, psi 

0 angle between normal line r and z-axis, degrees or radians 

\|r angle between normal line p and y-axis, degrees or radians 

Subscripts: 

i,j,k grid system 

x,y,z coordinate axes 

v, w,0 bay displacements 

o initial value 

n cycle of iteration 

A prime refers to the principal shear axes and 2 primes refer to the 
principal inertia axes. A bar over a symbol indicates an average value 
at the center of a bay. 
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Notation 

The notation employed is illustrated in figure 6. The system 
adopted for designating parts of the structure is as follows: 

Bulkheads divide the length of the structure into a number of bays. 
The subscript i is used to designate a given bulkhead or the bay 
between the i - 1 and ith bulkheads . 

The stringers and panels in a given cross section form the basis 
of a grid work which can be used to designate these elements. These 
grid lines are not necessarily straight, parallel, or perpendicular but 
follow the panels. Those grid lines that are approximately parallel to 
the z-axis are designated by the subscript jj those approximately 
parallel to the y-axis by the subscript k. 

With this system, points and stringers can be uniquely located as 
follows: 

The point on the ith bulkhead at the intersection of the jth and kth 
grid lines is .designated by the subscripts i,j,k. 

The stringer in the ith bay at the intersection, of the jth and kth 
grid lines is designated by the subscripts i,j,k. 

In order to locate a panel, the grid line on which it lies must be 
known. This notation consists of underlining the appropriate subscript} 
for example: 

The panel in the ith bay on the jth grid line and between the k. - 1 
and kth grid lines is designated by the subscripts i,j_,k. 

The panel in the ith bay on the kth grid line and between the j - 1 
and jth grid lines is designated by the subscripts i,j,k. 

The grid lines and bulkheads axe numbered such that the numbers 
increase in the directions of. the positive coordinate axes. 


Stress-Strain Relationships 

The shear strain in a given panel is constant and is defined by 
the following relationships which depend upon the location of the panel: 
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7 i>j,k 



= 2bTT( Ul ^> k + Ui - 1 .J. k " u i,J-l> k " »i-l,3-l,k) + 

J 

Av i Aw^ A0j_ 

LT cos + ET sln ^J< k ' ET r J< k 


(Ala) 


7 i,^ k ” 



2hJ~^( u i,3,k + u i-l,j,k - u i,j,k-l ’ u i-l,J,k-l) 


Avi Awi 

LT sln * ET C0S 


, A0i 

*J,k + lT Pj,k 


(Alb) 


When the shear strains are being computed, the normal distances r and 
p must be given their proper signs. 


The constant shear stress produces a linearly varying strain in 
the stringer and its average value at the center of the bay is 


_ u i,J,k “ u l-l,J,k 
Li 



(A2) 


Note that the thermal expansion is included in the relationship between 
stringer stress and strain. 


Equilibrium of Individual Stringers 

If a half -bay length of stringer on each side of point (i,J,k) is 
isolated, the force system of figure 7 is obtained, and the following 
equilibrium equation can be written: 

/ rtlA + /rtL\ _ / nrtL \ /rtL\ _ /jtL\ 

V 2 /i, j,k \ 2 /i,j+l,k \ 2 /i ,^,k \ 2 /i,^,k+l \ 2 /i+l,j,k 

P i,j,k + ^)i,j,k = 0 


(A3) 
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Substituting equations (Al) said (A2) into equation (A3) yields the 
following equation for the stringer displacement of point (i,j,k) in 
terms of the displacements of the adjacent bays and stringers: 
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where 



Equation (A4) involves no assumptions regarding equality of struc- 
tural dimensions, temperatures, or elastic properties about point (i,j,k) 
If any element is missing, the associated stiffness goes to zero and the 
general equation is still applicable. Since AE and Gt always appear 
as products, the variation of elastic properties with temperature is 
equivalent to changes in the stringer areas and sheet thicknesses of the 
effective structure. Furthermore, the thermal -expans ion terms appear in 
the same manner as axial loads applied to the stringers. Thus, if 
desired, the effects of a nonuniform temperature distribution can be 
determined by applying a set of equivalent loads to a new effective 
structure . 


Bay Shear and Torque Equilibrium 

The equations for the bay displacements (v,w ,0) can be obtained 
from equilibrium of the shear forces on the bay. cross section 


<Vi - EE 

j k 


(rtb cos 0)i,j,k ~ ( Tth sin = 0 


(*?>i ‘ Y Y Sin + (Tth C ° S = 0 


j k 


(A5a) 


(A5b) 



(A5c) 


Substitution of equations (Al) for the shear stresses in equations (A5) 
results in 
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Avi + (jyw'j i Awi - (^Qvj i A0j. - ^ + 

)i,j+l,k ( 


ILI^KjA 

3 k 


Gt cos 0 
2 


/Gt cos 0 
2 


i,J,k 


/Gt sin ilr \ 

1 2 i,i, 


k+l 


Gt sin t 




|jw) i Avi + ^Jww^ i Awi - ^J0v) i A0i = ^F z jj. + 

V~V~ / \P/Gt sin 0\ /Gt sin 01 

)_)_Kj,k + « 1 - 1 ,j,k) ( 2 A, J+1>k - ( 2 A,j,i 


J k 
/Gt cos 


/Gt cos ^ \ 

V 2 /i,J, 


k+l 


/Gt cos 


\ 2 

/ 




(jev)i ati - (j8w)i awi + (jee)i ASi = (m x )i - 

(Otr) 


i,J,k + u i-l,j,l 

j k 

/Gt£\ + /Gtp\ 

' 2 ' 2 '1,4,1 


l 1 2 • (n.« 


where 



/Gth sln 2 t 
l L 




*) -f 


Gth sin ^ cos ^ 


(A6a) 


(A6h) 


(A6c) 
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^Gth cos% j 




/Gthp sin ^ 



. Wi - , * • 

j k L ' J '- 


Gthp cos ^ 
L 






Equations (A6) can be simplified by eliminating the coupling terms 
if the axes used in the computations are the principal shear axes of the 
cross section. These axes are defined such that 


J vw’ 


= <v 


= «W 


= o 


(A7) 


The relationship between the location and orientation of points and 
panels in two systems of coordinates, arbitrary axes (x,y,z) and the 
principal shear axes (x'^y'^z 1 ), is shown in figure 8 and given by the 
following equations: 


- m')sin X' + (y - Z')cos X' 

(A8a) 

- m' )cos X' - (y - Z')sin X' 

(A8b) 

II 

(A9a) 

II 

-er 

(A9b) 

= r + Z' sin 0 - m' cos 0 

(AlOa) 

= p — Z 1 cos ^ - m' sin ♦ 

(AlOb) 


Then the location of the principal shear axes is 


tan 2X' 


2J. 


vw 


Jw 


J ww 


(Alla) 
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^vv^ww “ ^vw 

(Allb) 


, “ ^vw*^0w 

m 2 
JyvJww ~ Jw 

(Allc) 

with respect to these axes. 


Jyv 1 - Jw 

cos 2 X' + J w sin 2 X' + 2J W sin X' cos X' 

(A12a) 

^ww ' = 

o 2 

sin X' + J w cos X' - ZJyy sin X' cos X' 

(A12b) 


J ee' = J ee + i,J ew m ' J ev 

(A12c) 


Fy' = F z sin X' + Fy cos X' 

(A13a) 


F z ' = F z cos X' - Fy sin X' 

(A13b) 


V «M* + m'Fy - Z'F z 

(A13c) 


When referred to the principal shear axes, the equations for the 
hay displacements become 


Aw i ' 


(jyv 


L • 

HH( u i,J,k 

J k 

+ u i-i,j,k) 

/Gt cos 0' 

V 2 

/Gt 

COS 

1L\ 

/Gt sin 

, /Gt sin P\ 

r 

2 


' \ 2 J 

i^,k+l V 

2 k 

l 1 

yww 


r 

( F °')i * 

V 

OI( u i,J,k 

J k 

+ U i-l,J,k) 

/Gt sin 0' 
1 2 

/Gt 

sin 

o 

t) 

/Gt cos \ 

/Gt cos 

\ 

d 


' 2 h 

L, J.,k+1 ' 

' %. 


i,J+l,k 


I 


(Al4b) 
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11 = (s^Mi ' ZD U1 - J > k + - 

\ J k 

for' \ _ fop' \ fop' \ 

\ 2 'i,J,k V 2 fo^k+l ' 2 'i,±A 


(Al4c) 


Bay Thrust and Moment Equilibrium 

The equations obtained from equations (A4) and (Al4) are sufficient 
in themselves to define completely the displacements of the structure. 
However, if the equations are solved by iteration, it is helpful to 
employ a periodic correction cycle based on the gross equilibrium of 
axial loads on the cross section 


( F *)i - EIE(5)i,j, k - o , 

(A15a) 

(^i - = o 

(A15b) 

Ni + f ■ 0 

(A15c) 


It can be shown that equations (A15) are satisfied by the solution 
of equations (A4) and (Al4)j however, they are not -likely to be satis- 
fied by the displacement values obtained from any given cycle of itera- 
tion. In reference 2 it was demonstrated that convergence of the 
iterative process can be speeded if the displacement values are periodi- 
cally corrected so that the stringer displacements satisfy 
equations (A15). 

The corrections applied to the stringer displacements' are a planar 
distribution over the cross section and are determined as follows: 

( u i,j,k)n+l = ( u i, j,k)n + Au i,J,k (Al6) 


where 


Au i,j,k = 8 «i + Pz,iyj,k + Py,i 2 j,k 
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Substituting equation (Al6) into (A15) yields 

(AE) i bu i + (EQ z ) i 3 z ^ i + ( E %) i Py,^ i ' = (LF x )i + 


hLE 

J k 


( AEaT )i,J,k - (^) 


i > J >k 


(AlTa) 


(EQyJ^Ui + ( E Iyz)j3 z# i + ( EI yy)i3y^i = (E%)i 


QZ (AEaTz )i,j,k - ( u i,j,k - u i-i,j,k)(^-) 

j k ' ' 'i 


y J>kJ 


(AlTb) 


(EQ z )6u^ + (EI Z z)i3 z ^i + (Ely Z )^3y^i = -(LM Z )^ + 

Li 5ZH (AE a Ty)i,j,k - (ui,j,k - Ui-i,j, k )(^) 


■ J k 


J>k_ 


(A17c) 


where 


(AE)i = }_I_(AE)i,j, k 

(Elyy)i = ^ } (AEz^)i 1 k 

J k 

j k 

(KQ Z ) 1 . L)_(AEy)i,J,k 
J k 

(EI zz )i = ^^(AEy 2 )^^ 
J k 

J k 

(Eiyz)i = YI7~ (AEyz) i,j;k 

J k 


These equations can be simplified by elimination of the coupling 
terms if the computations are referred to the equivalent principal 
inertia axes of the cross section. These axes are referred to as equiva- 
lent because the variation of modulus of elasticity over the cross sec- 
tion is taken into consideration. These axes (x",y",z") are defined 
such that 


= 0 


EQ Z " = EQy” = EIy Z " 


(A18) 
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and then the following relationships are applicable: 
y" = (z - m" )sin X" + (y - Z")cos X 
z" = (z - m")cos X" - (y - Z")sin X 


tan 2X" = 


2(EI yz - AEZ"m") 

(Elyy -AEm" 2 ) - (ei zz - AEZ" 2 ) 



EQ Z 

AE 


(A19a) 

(A19b) 

(A20a) 

(A20b) 



(A20c) 


Elyy" = (Elyy - AEm."' 2 ) c o s : 2 '' + (ei zz - AEZ" 2 )sin 2 X" - 

. 2 (Elyz - AEZ"m")sin X" cos X" (A21a) 

EI ZZ " = (Elyy - AEm" 2 )sin 2 X" + (ei zz - AEZ" 2 ) cos 2 X" + 

2(EIy2 - AEZ"m")sin X " cos X" (A21b) 

F x " = F x * (A22a) 

My" = M z sin X" + My cos X" (A22b) 

Mz" = M z cos X" - My sin X" (A22c) 


A further simplification of the correction cycle can be made by 
eliminating the load and temperature terms on the right-hand side of 
equations (A17)- This elimination can be accomplished by iterating the 
difference between the exact solution and one which satisfies statics 
(equations (A15)) but not necessarily continuity. The iteration of dif- 
ferences has an additional advantage in that smaller numbers, and con- 
sequently less work, are required to obtain a solution. 
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An examination of equations (A17) indicates that they can be 
satisfied by a planar distribution of strain corresponding to the 
elementary^ analy sis of reference 1. Then the initial values cf stringer 
displacements u can be defined as follows: 

u i,j,k = u i-l,J,k + ^o^i^k (A&3) 

where 

( e oL)i,j,k: * (8Ui") 0 + (P z ,i”)'oyj,k" * <IV,l'Vj,k" 

and, with respect to the equivalent principal inertia axes, 

(**")o = (^).|C F x ,, )i +IZ^(AEaT) lrJA (A2ka) 

Oz,i”)o = ^7) i ^ M 'z ,, )i + [LWliJ, 

(3y,i”)o = + ^~^~^ AEaTZ '^ 1, - 1,k (A24c) 

The corresponding values of the bay displacements are obtained from 
equations (Al4) and (A23). 

Then the correction- cycle equations applicable to the iterated 
differences are as follows: 

8u i" = -(m) ± 5ZZj Au i,j,k - Au i-1 ,j,k)( AE )f,j,k (A25a) 

j k 

P z,i" = '(eTt) LD Au i,J,k - ^i-l,J,k)( AE y")i,j,k (A25b) 
' zz i j k 

P y,i" = "(eI^). - ^i-l,j,k)(AEz")i,j,k <A25c) 
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Equations (A25) provide corrections to the stringer displace- 
ments u only. These corrections remove any unbalanced moment or 
thrust on the cross section but add unbalanced shear forces which are 
removed by correcting the bay displacements ( v,w,0) . The corrected bay 
displacements are obtained from the corrected stringer displacements 
by application of equation (Al4). These two operations constitute the 
complete correction cycle that brings the stringer loads into equilibrium 
with the external loads without changing the shear stress in any panel. 
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TABLE VII.- CORRECTION CYCLE 


(a) Load problem. 
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Figure 1.- Typical multicell, stiff ened-shell wing structure. 



h — 3 at 10 = 30 


Figure 2.- Idealized -structure' used in illustrative example. (Stringer 
areas and skin thicknesses are listed in tables I and II.) 
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(a) Spanwise distribution of upper-surface stringer . stress . 



Figure 5-- Calculated stress distribution in the idealized structure. 
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Figure 5 • - Continued . 
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Figure 5*- Concluded. 







